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Abstract 

We present estimates of the masses of light quarks using lattice data. Our main 
results are based on a global analysis of all the published data for Wilson, Sheikholeslami- 
Wohlert (clover), and staggered fermions, both in the quenched approximation and 
with rif — 2 dynamical flavors. We find that the values of masses with the various 
formulations agree after extrapolation to the continuum limit for the rif — theory. 
Our best estimates, in the MS scheme at fi = 2 GeV, are m = 3.4 ± 0.4 ± 0.3 MeV 
and m s = 100 ± 21 ± 10 MeV in the quenched approximation. The rif = 2 results, 
m = 2.7 ± 0.3 ± 0.3 MeV and m s = 68 ± 12 ± 7 MeV, are preliminary. (A linear 
extrapolation in n; would further reduce these estimates for the physical case of three 
dynamical flavors.) These estimates are smaller than phenomenological estimates based 
on sum rules, but maintain the ratios predicted by chiral perturbation theory (%PT). 
The new results have a significant impact on the extraction of e'/e from the Standard 
Model. Using the same lattice data we estimate the quark condensate using the Gell- 
Mann-Oakes-Renner relation. Again the three formulations give consistent results after 
extrapolation to a = 0, and the value turns out to be correspondingly larger, roughly 
preserving m s (ipip). 
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1 INTRODUCTION 



The masses of light quarks m u , m d , and m s are three of the least well known parameters of the 
Standard Model. The range of values listed in the particle data book [1[], 2 < m u < 8, 5 < 
rrid < 15, and 100 < m s < 300 (evaluated in the MS scheme at ll = 1 GeV), is indicative of the 
large uncertainty. These quark masses have to be inferred from the masses of low lying hadrons. 
Theoretically, the best defined procedure is chiral perturbation theory (xPT) which relates the 
masses of pseudoscalar mesons to m u , rrid, and m s . Recall that the lowest order chiral Lagrangian 
is 

f 



C = — tr 



{dpEduTt) + 2^(ME + ME 1 ") (1) 



The presence of the overall unknown scale /i implies that only ratios of quark masses can be 
determined using xPT. The predictions from xPT for the two independent ratios are [0. [|. [| 



Lowest order Next order 
2m s /(m u + m d ) 25 24.4(1.5) 

m u /m d 0.55 0.553(43) . 



One notes that the change from lowest to next order is insignificant. Unfortunately, the prospects 
for further improvement within xPT are not good. To get absolute numbers, one has to combine 
the xPT analysis with sum rule or with some other model calculations. The two most recent and 
up-to-date versions of the sum rule analysis has been performed by Bijnens, Prades, and de Rafael 
giving m u + rrid = 12(2.5) MeV and by Jamin and Munz |J giving m s = 178(18) MeV. It is 
typical of sum rule analysis to quote results at ix = 1 GeV, while lattice results are presented at 
ix = 2 GeV. To facilitate comparison, we translate the sum-rule values to lm — 2 GeV, whereby 
m u + m d = 9.4(1.8) MeV and m s = 126(13) MeV using A^L = 300 and 380 MeV respectively. 
These values are about one a above the quenched lattice results we present below. However, as 
discussed in J7|, a reanalysis of the sum- rules calculations show that the uncertainties are large 
enough to preclude any serious disagreement with the lattice results. To improve upon the sum- 
rule estimates requires hard to get at experimental information on the hadronic spectral function. 
Thus, we believe that lattice QCD offers the best approach to determining these quantities, and 
this paper presents an analysis of the current data. 

In lattice QCD the quark masses are input parameters in the simulations. Their values are 
determined by tuning the masses of an equal number of hadrons to their physical values. One 
additional hadron mass is needed to fix the lattice scale, or equivalently the value of a s . In our lattice 
simulations we do not include electromagnetic effects, so we can only calculate the isospin symmetric 
light quark mass m = (m u + rrid)/2. The hadrons we use to fix the lattice scale are M p = 770 MeV 
(one could instead, for example, use f n = 131 MeV, but its determination is less reliable and has not 
been reported for all the data we use), while to fix m and m s we study the behavior of pseudoscalar 
and vector meson masses as a function of the quark mass. (For experimental numbers we use the 
isospin averaged values M n = 137 MeV, M K = 495 MeV, M = 1020 MeV and M* K = 894 MeV.) 
We have chosen to use these pseudoscalar and vector mesons since the corresponding correlation 
functions have been measured with the smallest statistical errors. Thus, in principal, quark masses 
can be determined without any ambiguity from these calculations. 

Current lattice simulations make a number of approximations because of which estimates have 
systematic errors in addition to statistical errors. Some of these issues have been discussed in the 
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reviews by Ukawa || and Gupta |J. In this paper we present an analysis of the cumulative world 
data as available in November 1996 for staggered, Wilson, and Sheikholeslami-Wohlert (c$w — 1) 
fermions, and both for the quenched and n/ = 2 theories. (The parameter csw is the strength of 
the clover term that is added to the Wilson action to remove the 0(a) discretization errors [H], |H|.) 
We present evidence showing that the largest uncertainty in the results arises from discretization 
(0(a)) errors and from the dependence on the number of dynamical flavors. We find that the 
discretization errors show the expected behavior, 0(a) for Wilson and 0(a 2 ) for staggered. For 
clover action one expects corrections of 0(a s a), however the current data show behavior similar 
to that with Wilson fermions. By combining the results from these three formulations we show 
that the discretization errors can be controlled. Thus the main remaining uncertainty is due to 
extrapolation in n/, because of which we are only able to extract rough estimates of the quark 
masses for the physical case of rif = 3. These are significantly smaller than the phenomenological 
estimates mentioned above, and lead to a large enhancement in the standard model prediction of 
e'/e. 

Lastly, we also present an analysis of the quark condensate, (ipip) , using the data for pseudoscalar 
meson mass and the Gell-Mann-Oakes-Renner relation. Again we find that the three lattice dis- 
cretizations of the Dirac action give consistent results after extrapolation to the continuum limit. 
The condensate turns out to be roughly a factor of two larger than phenomenological estimates 
such that the renormalization group invariant quantity m{%p%p) is preserved. 

This paper is organized as follows. In Section ^| we summarize the definition of the quark mass, 
the relation between the lattice and continuum results, and the 2-loop running. The reorganization 
of the lattice perturbation theory a la Lepage and Mackenzie is discussed in Section |3|. A brief 
description of the lattice data used in this analysis is given in Section f|. Analyses of m and m s are 
presented in Sections |5| and We discuss the variation of quark masses with the clover coefficient 
in Sheikholeslami-Wohlert action in Section [7[ A general discussion of the systematic errors is given 
in Section || In Section || we compare our results with earlier calculations and discuss the impact 
on e'/ e in Section [10]. The calculation of the quark condensate is presented in Sections [□]. We end 



with our conclusions and future outlook in Section O 



2 DEFINITION OF QUARK MASSES 



There are two ways of calculating the MS mass for light quarks at scale /x from lattice estimates. 
The first is 

m Ms(/ i ) = Z m (fia)m L (a) 

1 - \ [8log(fia) - C m ]|m L (a) 

1 - \[8log(fia) - (C m - tad)} \ {m L (a)X} (2) 



where Z m = l/Z$ is the mass renormalization constant relating the lattice and the continuum 
regularization schemes at scale /i, and A = g 2 /16Tc 2 . The 1-loop perturbative expressions for Z m are 
given in Table |l| For Wilson (or Sheikholeslami-Wohlert (SW) clover) type of fermions, the lattice 
estimate of the quark mass, defined at scale q*, is taken to be m^(g*)a = (1/2k — 1/2k c ). We have 
also carried out the analysis using the alternate definition Ln(l + (1/2/t — 1/2k c )). The change in 
individual estimates is < 3% for data at a < 0.5 GeV -1 and negligible in the extrapolated values. 
Since the analysis with this 0(a) improved definition is much less transparent we use the simpler 



2 



form. For staggered fermions rriL(q*) = m Q , the input mass. The last form in Eq. ||] shows the result 
of applying tadpole improvement as discussed in Section [| 

The second method is to use the Ward Identity for the renormalized axial vector current, 
O^ZaA^ = (mi + m2)ZpP + 0(a), where mi and are the masses of the two fields in the 
bilinear operators. The quark mass, to 0(a), is then given by the following ratio of correlation 
functions 

(roi+m , )= jjWW)), (3) 
V ; Z P (P(r)J(O)) y 1 

where J(0) is any interpolating field operator that produces pions at p = 0. Here, Za and Zp are 
the renormalization constants for the axial and pseudoscalar densities that match the lattice and 
continuum theories at scale \x = q*. The important point to note is that since the Ward identity 
is valid locally, the value of m should be independent of r. Violations of this relation are signals 
for 0(a) errors in the currents. This has been exploited by the Alpha collaboration to remove the 



0(a) discretization errors non-perturbatively in the clover action and currents [|Tl| . 

At long time separation (r large) one can assume that the 2-point correlation function is de- 
scribed by the asymptotic form ~ e~ m " nT . One can then write Eq. [3] as 

(m | m ) Z _A m (Mr)J(O)} ( ) 

{mi + m2 >- Zp™* (P(r)J(O))- (4) 

There are very few calculations of m q using these two versions of the Ward Identity. Also, at this 
stage the perturbative and non-perturbative estimate of Zp do not agree ||12fl. For these reasons we 



postpone the discussion of this method to a later work [|13| . 

Once m-jjg has been calculated at the continuum scale /i, its value at any other scale Q is given 



by the two loop running [O] 



m(Q) fg 2 (Q)\^ ( x + g 2 (Q) -g'M^ift-ToA j | |r)) 



m(ji) \g 2 (n)J V 16tt 2 v 2/5, 



o 



where in the MS scheme for iV colors and rif flavors 

A) 



11N -2n f 



34 Ar2 10 Ar N 2 -l 

3 X = —N 2 Nn f n f , 

H 3 3 N 

a N 2 -l 

' 3 2N v 2N ' 3 2N K ' 

We will quote our final numbers at Q = 2 GeV as we do not feel confident using the 2-loop relations 
to run to smaller scales. 



3 Renormalization Constants 



Our final results are stated in the MS scheme. The required 1-loop renormalization constants for 
matching the lattice and continuum operators for Wilson, clover, and staggered versions of the 
lattice discretization are collected in Table [I] using results derived in [TJJ, [15], [16]. We reorganize 
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Table 1: Renormalization constants in the MS scheme before tadpole subtraction for staggered |15| ], 
and Wilson and clover (c$w — 1) []16| fermions. Here A = g 2 /167r 2 . 





Staggered 


Wilson 


Clover(csw = 1) 


Z A 

z P 

Z m 


1 

1/ Z m 

1 - A(81n(/ia) - 52.288) 


1 - 21.06A 
1 + A(81n(/ia) - 30.13) 
1 - A(81n(/xa) - 17.27) 


1 - 18.39A 
1 + A(81n(/xa) - 29.84) 
1 - A(81n(/ia) - 25.75) 



lattice perturbation theory using the Lepage-Mackenzie prescription |jr7|| . This prescription involves 
four parts: the renormalization of the quark field Z^, and the quark mass, the removal of tadpole 
contribution from the 1-loop operator renormalization, the choice of a s , and an estimate of the 
typical scale q* characterizing the lattice calculation. 



We use the following Lepage-Mackenzie definition of the strong coupling constant |T7|, [18], [19] 



1 47T 

-\n{-Trplaq) = — a v (3.41/a) (1 - (1.191 + 0.025n f )a v ) (Wilson), 

o o 

47T 

= —a v (3Al/a) (1 - (1.191 + 0.070n/)a y ) (Stag.), (7) 



from which the MS coupling at scale 3.41/a is given by 

2 



ajrs^Al/a) = av-(e 5/6 3.41/a)(l + -a v ). (8) 

7T 



The value of otj^-g at any other scale is then obtained by integrating the standard 2-loop /5-function 
for the appropriate number of flavors. The results are given in Tables ^|, BL and |6[. Ideally, 
GeV) should be independent of (3. From these tables one can see the variation with /3, 
and the extent to which the variation in 1/a at fixed (3 feeds into a^gg(2 GeV). Reorganizing the 
lattice perturbation theory in terms of ol-jj^ has the advantage that the continuum and lattice a s is 
the same when matching the theories at scale = q*. We call this matching procedure "horizontal" 



matching [ 20| . 



An estimate of q* is not straightforward since Z m is logarithmically divergent. So we appeal to 
the general philosophy of tadpole improvement, i.e. it is designed to remove the short-distance lattice 
artifacts. Once these have been removed, the typical scale q* of the lattice calculation becomes less 
ultraviolet. Thus, our preferred scheme is one in which q* = u = 1/a, which we call TAD1. To 
analyze the dependence of the results on q* we also investigate the choices q* = u = 0.5/a,2/a 
and re/a. A more detailed description of these schemes and our implementation of the matching 
between the lattice and continuum is given in . 



In all methods of calculating the quark mass, Eqs. 0, |3], ^, the normalization of the quark fields 
does not enter, or cancels between the different currents. Thus, one does not have to consider this 
factor. The quark mass gets scaled by the tadpole factor X, i.e. —>■ Xm^, where we use the 
non-perturbatively determined estimate for X. 



The tadpole factor X for Wilson and clover fermions is chosen to be 8k c pl| . It has the 
perturbative expansion X = 1 + tad A where 

tad = 17.14 , (Wilson) , 

tad = 10.66 , (Clover (c sw = 1)) • (9) 

For Wilson fermions at, say, j3 = 6.0, A = 0.0153 in the TAD1 prescription, so the perturbative 
value 8k c = 1.262 is almost identical to our non-perturbative estimate 8k c = 1.257 ||22|| . The 
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Lepage-Mackenzie reorganization, as shown in Eq. |2|, therefore seems benign, nevertheless, we in- 
clude tadpole subtraction as its purpose, in general, is to reorganize lattice perturbation theory to 
make the neglected higher order corrections small. The entries in Table |l] show that after tadpole 
subtraction the 1-loop corrections are indeed smaller for Wilson and clover fermions. However, note 



that the 1-loop results presented in |Tj|, and used here, show that the perturbative correction grows 
with c S w 

For staggered fermions we define the tadpole factor X as the inverse of the fourth root of the 
expectation value of the plaquette 

X = Uq 1 = plaquette~ 1/4 = (1 + 13.16A) . (10) 

The agreement, at /3 — 6.0, between the perturbative value 1.20 and the non-perturbative value 
1.14 is still reasonable in TAD1, even though the plaquette is designed to match at q* ~ 3.41/a as 
shown in Eq. [7]. The point to note is that for any reasonable choice of X, the correction C m — tad 
is still very large for staggered fermions . Thus one might doubt whether the 1-loop perturbative 
result for Z m is reliable. It turns out, as we show later, that using the expressions given in Table [1] 
give results that, in the a = limit, agree with those obtained using Wilson-like fermions. 

The second noteworthy outcome of our analysis is that the dependence of the results on the 
choice of q*, and whether or not one does tadpole subtraction, are small. We find that m increases 
with q*, and the most significant variation is in the quenched staggered data (~ 5%). The difference 
between tadpole subtraction and no tadpole subtraction is significant only at strong coupling (~ 3% 
at P = 6.0 for the quenched theory, and < 3% for f3 > 5.4 for rif = 2), and even smaller for different 
choices of tadpole factor. Based on the above estimates we assume that for a < 0.5 the uncertainty 
in relating the lattice results to those in a continuum scheme is under control at the 5% level. We 
consider this variation negligible since for many of the data points it is small compared to even the 
statistical errors. Furthermore, this uncertainty is one aspect of the discretization errors, and our 
estimate of the errors associated with the extrapolation to a — > incorporates it. To summarize, 
all of our lattice estimates are quoted using the TAD1 scheme, and the uncertainty associated with 
this choice is included in our quoted extrapolation error. 



4 LATTICE PARAMETERS OF DATA ANALYZED 



The list of calculations from which we have taken data are summarized in Tables |2| and 
quenched Wilson data are taken from Refs. [li 
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m 3i, p 



], quenched clover from 
2) staggered fermions from [19], 



The 

|28fl , quenched staggered from 
2) Wilson from @, |35 



dynamical 

|39|. All these calculations use the simple 



and dynamical (rif 

Wilson (plaquette) action for the gauge fields. Thus the difference between staggered and Wilson- 
like fermion data for fixed nf is a reflection of the difference in 0(a) errors between the various 
discretization schemes. 

In principle, if lattice simulations could be done with three dynamical flavors and with realistic 
parameters, one could tune the three quark masses to reproduce the spectrum and thereby determine 
their values. However, since this is not the case, the strategy we use is as follows. For each value of 
the lattice parameters (/3, rif, fermion action) we fit the data for M% and M p as a linear function 
of the quark mass 



Ml = A 7T + B n (m 1 + m 2 )/2 
M p = A p + B p (m 1 + m 2 )/2 
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Table 2: Lattice parameters of the quenched data used in the global analysis 
the type of fermion action (jy=Wilson, 5'=staggered, C=clover with csw = 
6/g 2 , the lattice size, the number of configurations in the statistical sample, 
masses in lattice units (k values) used in the fits for staggered (Wilson and 
lattice scale 1/a GeV extracted from M p , and the values of a ~jjg at g* = 1/a 



We list the reference, 
1), the coupling (3 = 
, the values of quark 
clover) fermions, the 
and 2 GeV. 



Label 



Ref. 



Lattice 



size 



#of 
Conf. 



Quark masses used in the fits 



Scale 1/a 
(GeV) 



a- 



l MS 
1/a) (2 GeV) 



Wl 

W2 

W3 

W4 

W5 

W6 

W7 

W8 

W9 

W10 

Wll 

W12 

W13 

W14 

W15 

W16 

W17 

W18 

W19 

W20 



H 

1 



25j 
1 

g 

m 
m 
n 

m 
m 
m 

n 



5.7 
5.7 
5.85 
5.85 
5.9 
5.93 
5.95 
6.0 
6.0 
6.0 
6.0 
6.1 
6.17 
6.2 
6.2 
6.2 
6.3 
6.3 
6.4 
6.4 



24 3 x 32 
24 3 x 32 
24 3 x 54 
16 3 x 32 
16 3 x 40 
24 3 x 36 
16 3 x 32 
32 3 x 64 
24 3 x 32 
24 3 x 54 
18 3 x 64 
24 3 x 64 
32 3 x 40 
24 3 x 48 
24 3 x 64 
24 3 x 64 
32 3 x 80 
24 3 x 32 
24 3 x 60 
24 3 x 64 



50 

92 
100 

90 

150 

210 

90 
170 

78 
200 
320 
100 
219 

18 
250 
110 
100 
128 

15 
400 



0.165, 0.167, 0.168 
0.165, 0.1663, 0.1675 
0.1585, 0.1595, 0.1605 
0.1585, 0.1600 
0.156, 0.157, 0.158, 0.1585 
0.156, 0.1573, 0.1581 
0.1554, 0.1567 
0.155, 0.1558, 0.1563 
0.155, 0.1558, 0.1563 
0.155, 0.1555, 0.1563 
0.153, 0.154, 0.155 
0.152, 0.153, 0.154, 0.1543 
0.1519, 0.1526, 0.1532 
0.1523, 0.1526, 0.1529 
0.1510, 0.1515, 0.1520, 0.1526 
0.1510, 0.1520, 0.1526 
0.150, 0.1505, 0.151, 0.1513 
0.1485, 0.1498, 0.1505 
0.1485, 0.1490, 0.1495 
0.1488, 0.1492, 0.1496, 0.1500 



1.431(27) 
1.422(24) 
1.958(114) 
1.741(76) 
1.987(48) 
2.000(37) 
1.941(58) 
2.338(43) 
2.204(70) 
2.423(146) 
2.154(64) 
2.629(59) 
2.755(48) 
2.735(172) 
2.914(89) 
2.934(121) 
3.260(88) 
3.092(68) 
3.628(416) 
4.095(163) 






246 





211 





246 





210 





213 





211 





213 





201 





205 





204 





201 





201 





198 





196 





192 





205 





192 





200 





192 





208 





192 





198 





181 





201 





176 





198 





173 





194 





173 





199 





173 





200 





166 





197 





166 





193 





159 





195 





159 





205 



SI 

S2 

S3 

S4 

S5 

S6 

S7 

S8 

S9 

S10 

Sll 

S12 

S13 

S14 

S15 



m 



m 
u 

M 

m 
n 



32j 



51 



5.7 
5.7 
5.85 
5.85 
5.93 
5.95 
6.0 
6.0 
6.0 
6.0 
6.2 
6.2 
6.4 
6.4 
6.5 
6.5 



24 3 x 32 
16 3 x 32 
16 3 x 32 
16 3 x 32 
20 3 x 40 
16 3 x 32 
24 3 x 32 
24 3 x 40 
32 3 x 64 
24 3 x 64 
32 3 x 48 
32 3 x 64 
32 3 x 48 
40 3 x 96 
32 3 x 64 
48 3 x 64 



50 
32 
160 
60 
50 
190 
60 
23 
200 
50 
23 
40 
24 
40 
100 
200 



0.02, 0.015, 0.01, 0.005 

0.02, 0.015, 0.01, 0.005 

0.025, 0.01 

0.01, 0.02, 0.03, 0.04 

0.01, 0.02, 0.03, 0.04 

0.025, 0.01 

0.04, 0.02, 0.01 

0.03, 0.025, 0.02, 0.015, 0.01 

0.01, 0.005, 0.0025 

0.02, 0.03, 0.04 

0.025, 0.015, 0.01, 0.005 

0.005, 0.01, 0.02 

0.015, 0.010, 0.005 

0.005, 0.01, 0.02 

0.01, 0.005, 0.0025 

0.01, 0.005, 0.0025 



0.951(78) 
1.312(75) 
1.340(31) 
1.571(33) 
1.645(25) 
1.843(78) 
1.911(84) 
1.917(34) 
1.855(41) 
2.569(78) 
2.616(91) 
3.467(376) 
3.429(78) 
3.962(127) 
3.811(59) 



0.246 
0.213 
0.213 
0.201 
0.198 
0.192 
0.192 
0.192 
0.192 
0.173 
0.173 
0.159 
0.159 
0.152 
0.152 



0.180 
0.180 
0.182 
0.183 
0.184 
0.187 
0.189 
0.189 
0.187 
0.189 
0.191 
0.192 
0.191 
0.192 
0.189 



CI 
C2 
C3 
C4 



m 
m 



6.0 
6.2 
6.2 
6.4 



18 3 x 64 
24 3 x 64 
18 3 x 64 
24 3 x 64 



200 
250 
200 
400 



0.1425, 0.1432, 0.1440 
0.14144, 0.14184, 0.14224, 0.14264 
0.14144, 0.14190, 0.14244 
0.1400, 0.1403, 0.1406, 0.1409 



1.867(81) 
2.604(131) 
3.082(384) 
3.951(182) 



0.192 
0.173 
0.173 
0.159 



0.187 
0.190 
0.204 
0.202 
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Table 3: Lattice parameters of the rif 
The rest is same as in Table [|. 



2 dynamical simulation data we use in the global analysis. 



Label 


Ref. 


/3 


Lattice 
size 


#of 
Conf. 


Quark masses used in the fits 


Scale 1/a 
(GeV) 


(1/a) (2 GeV) 


W 2 1 
W 2 2 
W 2 3 
W 2 4 
W 2 5 
W 2 6 


m 

[34| 
[36| 

m 
\m 


5.3 
5.4 
5.5 
5.5 
5.6 
5.6 


16 3 x 32 
16 3 x 32 
16 3 x 32 
16 3 x 32 
16 3 x 32 
16 3 x 32 


417-484 

14- 15 

15- 27 
400-669 

32-45 
100 


0.167, 0.1675 
0.160, 0.161, 0.162 
0.158, 0.159, 0.160 
0.1596,0.1600, 0.1604 
0.156, 0.157 
0.156, 0.157, 0.1575 


1.890(110) 
1.467(113) 
1.812(97) 
2.179(84) 
2.332(209) 
2.379(77) 


0.266 0.259 
0.254 0.223 
0.227 0.219 
0.226 0.234 
0.213 0.226 
0.213 0.227 


5 2 1 
5 2 2 
5 2 3 


[39| 




5.6 
5.7 
5.7 


16 3 x 32 
16 3 x 32 
20 3 x 20 


200-400 

-50 
150-160 


0.025, 0.01 

0.025, 0.02, 0.015, 0.01 
0.02, 0.01 


1.798(73) 
2.341(78) 
2.236(99) 


0.232 0.223 
0.216 0.230 
0.216 0.226 



where mi and m 2 are the masses of the two quarks. (We use M n as shorthand for pseudoscalars, 
and M p for vector mesons.) These fits are made assuming that the data at each value of the quark 
mass are independent, and the statistical error estimates quoted by the authors are used in the 
X 2 minimization procedure. Note that for Wilson-like fermions we make fits to 1/2k whereby 
defines k c and A p is the intercept at 1/2k = 0. For staggered fermions A n should be zero due to 
chiral symmetry. However, in our analysis we leave it as a free parameter and use it to define the 
zero of the quark mass. This leaves B n , A p , B p from which we determine the three quantities that 
we are interested in; the scale using M p , m using M^/M^, and m s in three different ways using 
M K , M K *, M$. Throughout the analysis we assume that is a pure ss state. 

The error analysis of the global data is somewhat limited since in most cases we only have access 
to the final published numbers. For this reason we cannot include two types of correlations in the 
data when making fits using Eq. [11]. Firstly, those between pseudoscalar and vector mesons at a 
given quark mass, and secondly those between the meson masses at different quark masses calculated 
on the same background gauge configurations as is the case for the quenched simulations. Thus, 
when making chiral fits we simply assume that the errors in the lattice measurements of meson 
masses are uncorrelated. Thereafter, the errors in A n , B n , A p , B p are propagated self-consistently to 
the final estimates. By comparing results in a few cases where a full error analysis has been done, we 
find that these shortcomings change estimates of individual points by less than even the statistical 
errors. The one exception is the results from QCDPAX collaboration [26 where the neglected 



correlations have a large effect. They quote significantly different values for A n , B n , A p , B p based 
on correlated fits. Using their correlated fit parameters give quark masses consistent with other 
estimates at the same couplings. Thus, we believe that our estimates based on a reanalysis of the 
global data are reliable. 

The input parameters and results ( a ]\jg> lattice scale 1/a, A n , B n , A p , B p , and the quark masses 
in lattice units and in MeV in MS scheme at 2 GeV) for each data set are given in Tables 0, [| |6], |, 
and [7|. (In the case of data from Ref. |30| , the analysis started with the results for A n , B n , A p , B p 
because the raw data have not been published.) We note that the various data points are obtained 
on lattices of different physical volumes, the statistical sample size varies significantly (the sample 
size is rather small in some cases), and the strategy for extracting masses varies from group to 
group. However, it turns out that in the calculation of quark masses, there is a cancellation of 
errors that make these differences much less significant than for example in the estimates of meson 
masses themselves. We substantiate some of these remarks by presenting a detailed analysis of 
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some of the systematic errors in Section || using data from our high statistics calculation presented 
in Ref. 0. 

Once we have the quark masses at different lattice scales, we analyze the dependence on the 
lattice spacing a by comparing Wilson, clover, and staggered results. Where data permit, we omit 
points at the stronger couplings (larger a) for the following two reasons. First, we use only the 
leading correction in the extrapolation to a = 0, and secondly, the perturbative matching becomes 
less reliable as (3 is decreased. The end result is that we find that the leading corrections give a good 
fit to the data, and in the a = limit the different fermion formulations give consistent results. 

To investigate the n/ dependence, we compare the rif = and rif = 2 results. We then assume 
a linear behavior in rif to extrapolate these two points to estimate the desired physical value. As 
we make clear in our discussions, the extrapolation in rif turns out to be the weakest part of our 
analysis. 

5 LIGHT QUARK MASS m = (m u + m d ) /2 

In order to extract light quark masses from lattice simulations we have used the simplest ansatz, 
Eq. [II], for the chiral behavior of hadron masses. The reason for this truncation is that in most 
cases the data for M n and M p exist at only 2 — 4 values of "light" quark masses in the range 
0.3m s — 2m s . In this restricted range of quark masses the existing data do not show any significant 
deviation from linearity. Using linear fits to M 2 data over a limited range in m q means that we 
can predict only one independent quark mass from the pseudoscalar data, and we cannot test the 
difference between tree level and NLO results predicted by analyses. The mass we prefer to 
extract using the pseudoscalar spectrum, barring the complications of quenched ;\TT which are 
expected to become significant only for m q < 0.3m s || fjOfl , is m because to extract m s using m K 
needs first an extrapolation to m in the light quark. While this choice avoids the question whether 
lowest order %PT is valid up to m s , one must bear in mind that the fits are made to data in the 
"heavier" range 0.3m s < m q ~ 2m s and then extrapolated to m. The bottom line is that to get m 
we linearly extrapolate the ratio M 2 /M 2 in m to its physical value 0.03166. 

The global quenched data are shown in Fig. [IJ and listed in Table f| for the tadpole subtraction 
scheme TAD1. The expected leading term in the discretization errors is 0(a) for Wilson, 0(g 2 a) 
for clover, and 0(a 2 ) for staggered. The fits, keeping only the leading dependence (for clover data 
we assume a linear dependence in a), give 

m(a) = 3.33(22) MeV [1 + 1.3(2) GeV a] X 2 /dof = 0.8 Wilson 

m(a) = 3.13(42) MeV [1 + 1.4(4) GeV a) X 2 /dof = 0.6 Clover 

m(a) = 3.27(04) MeV [1 + (0.20(11) GeV a) 2 } x*/dof = 4.6 Staggered . 

For Wilson and staggered fermions we fit to (3 > 5.93 (a < 0.5 GeV -1 ) to minimize the effect of 
higher order corrections. The staggered data show a small rise for (3 > 6.0 after the initial fall. With 
present data it is not clear whether this is due to poor statistics, the decreasing physical volume 
of the lattices used, or higher order terms in a and g 2 . As a result, the fit with only the lowest 
order correction is not able to capture this trend in the data. This is reflected in the high \ 2 jdoj . 
Thus we consider, as an alternate estimate of the continuum result, the mean, m = 3.4(1) MeV, 
of the estimates at the two highest (3. A second alternative is to fit the data for (3 > 6.0. This 
give m = 3.53(6) MeV with a y^jdoj = 0.95. The corresponding fit for Wilson fermions gives 
m = 3.45(23) MeV with a X*/dof = 0.9. 

Given the size of 0(a) corrections in the case of Wilson fermions, we have also fit the data 
including a quadratic correction as shown in Fig. |2|. Unfortunately, we find that even the sign of 
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Table 4: Results for the fit parameters k c for Wilson-like and A w for staggered fermions, B w , A p 
and Bp, defined in Eq. [TlJ and the quark masses for the quenched simulations. Both the lattice 
values, ma and m s (M<^)a, as well as the physical values, m(MS,2 GeV) and m s (MS,2 GeV, M$) 
in MeV, are given. The lattice parameters and the references to the original work are given in 
Table ^| Chiral fits for labels marked by a * have x 2 /d°f — 3.0. 
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Table 5: The results for the rif — 2 dynamical simulations listed in Table |3|. The notation is same 
as in Table ||. 



Label 



Kc/ A n 



B n 



m 

ma x 10 5 



m 



m s (M^) 
m°a x 10 3 m, 



W 2 1 
W 2 2 
W 2 3 
W 2 4 
W 2 5 
W 2 6* 



.167928(07) 

164866(288) 

.161335(64) 

.161218(08) 

.158393(47) 

.158506(13) 



12.471(365) 
6.287(440) 
4.979(225) 
4.499(134) 
4.426(285) 
3.903(87) 



-41.186(4356) 
-10.566(1450) 
-11.542(1239) 
-14.173(1801) 
-14.249(2216) 
-13.211(823) 



13.967(1455) 
3.655(466) 
3.860(393) 
4.683(576) 
4.617(694) 
4.290(258) 



42(05) 1 

139(20) 2 

115(11) 2. 

88(07) 2 

78(13) 2 

85(05) 2 



06(06) 
53(20) 
64(15) 
52(11) 
37(21) 
65(09) 



10(02) 
48(09) 
37(06) 
25(04) 
24(06) 
25(02) 



25(3) 
88(11) 
85(9) 
73(9) 
73(11) 
79(5) 



<S 2 1* 

S 2 2 

<S 2 3 



-0.00025(1) 
0.00012(7) 
-0.00078(11) 



6.951(66) 
6.041(98) 
5.573(162) 



0.423(17) 
0.323(11) 
0.345(14) 



8.667(777) 
9.120(501) 
7.320(762) 



84(07) 2 
57(04) 2 
67(06) 2 



93(12) 
62(09) 
95(14) 



17(02) 
12(01) 
15(02) 



58(5) 
57(3) 
66(8) 



the quadratic term changes depending on whether we include all the points with (3 > 5.93 or only 
those with (3 > 6.0. Also, one cannot distinguish between the three fits on the basis of x 2 /dof 
which lies between 0.8 — 0.96. The lack of stability of the quadratic fits indicates that possible 
0(a 2 ) corrections cannot be determined with the current data. For this reason we use the results 
obtained using the lowest order fit as our best estimate. 

The clover (c$w — 1) fermion results are surprisingly close to the Wilson values even though 
the discretization errors should be smaller, i.e. 0(g 2 a). In particular the variation with f3 is similar 
to that for Wilson fermions, however it should be noted that this is based on data at only three 
values of /3. At present, looking at the trend in the data, the best we can do is to assume a linear 
behavior. 

The bottom line is that the results, in the a = limit, from all three formulations turn out to 
be in surprisingly good agreement with each other. Thus, for our best estimate of the quenched 
value we take an average of the f3 > 6.0 numbers 

m(MS, 2 GeV) = 3.4 ± 0.4 ± 0.3 MeV (quenched). (12) 

where the first error estimate is the largest of the extrapolation errors and covers the spread in the 
data. The second error is a 10% uncertainty due to the lattice scale as discussed in Section [| 

To analyze the rif = 2 dynamical configurations, we have restricted ourselves to data with 
^valence = ^sea- The main limitation we face is that the data have been obtained at very few values 
of lattice spacing, and the statistics and lattice volumes are smaller than in quenched simulations. 
The pattern of 0(a) corrections in the present unquenched data, shown in Fig. |3| and Table |5|, is not 
clear. In fact, as discussed in Section |8], we cannot even ascertain whether the convergence in a is 
analogous to the quenched case, i.e. from above for both Wilson and staggered formulations. The 
strongest statement we can make is qualitative; at any given value of the lattice spacing, the rif = 2 
data lies below the quenched result. Taking the existing data at face value, we find that the average 
of the Wilson and staggered values are the same for the choices (3 > 5.4, (3 > 5.5, or (3 > 5.6. Since 
for /? > 5.5 there are two independent measurements that agree, our current estimate is taken to 
be the average of these data 

m(2 GeV) = 2.7 ± 0.3 ± 0.3 MeV (n f = 2 flavors) , (13) 
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where the first error estimate is the spread in the data. To obtain a value for the physical case of 
nf = 3, the best we can do is to assume a behavior linear in rif. In this case a linear extrapolation 
of the nf = and 2 data gives 

m(2 GeV) « 2.4 MeV (n/ = 3 flavors). (14) 

It is obvious that more lattice data are needed to resolve the behavior of the unquenched results. 
However, the surprise of this analysis is that both the quenched and n/ = 2 values are small and 
lie at the very bottom of the range predicted by phenomenological analyses ||. 



6 THE STRANGE QUARK MASS m s 

We determine m s using the three different mass-ratios, M^/M 2 , Mk*/M p , and M^/Mp. As men- 
tioned above, using a linear fit to the pseudo-scalar data constrains m s (M^)/m = 25.9. Unfor- 
tunately, the data are not good enough to include terms of order m 2 or the non-analytical chiral 
logarithms when making fits. Thus we cannot improve on the lowest order results 



m s (MS, n = 2 GeV, M K ) = 88(10) MeV (Quenched), 

m s (MS, /i = 2 GeV, M K ) = 70(8) MeV (n f = 2) , (15) 

when extracting m s from M\jM 2 . Using the vector mesons and M<f, gives independent esti- 
mates. To illustrate the difference, a comparison of m s (M^) and m s = 25.9m is shown in Fig. |] for 
the quenched Wilson theory. 

The quenched Wilson and staggered data, and the extrapolation to a = 0, for m s (M^) are shown 
in Fig. |5|. The Wilson data again show large 0(a) corrections and we make a linear fit to the data 
at (3 > 5.93. The staggered data converges from below, and anticipating 0(a 2 ) corrections we fit 
linearly in a 2 . The clover data are shown in Fig. |6|. The parameters of the fits are 

m s (M^) = 94(18) MeV[l + 1.9(7) GeV a] X 2 /dof = 0.4 Wilson 

m s (M ) = 62(31) MeV[l + 3.9(2.5) GeV a] x 2 /dof = 0.34 Clover (16) 

m s (M ) = 99(4) MeV[l - (0.48(15) GeV a) 2 ] x 2 /dof = 2.0 Staggered . 

Compared to m, Wilson and clover fits have much larger errors. This is due to the fact that 
the lattice measurements of the vector mass have much larger statistical errors compared to the 
pseudoscalars. Also, the spread between different groups is larger, reflecting the differences in the 
strategy to extract the masses from the 2-point correlation functions. The fit to the staggered data 
for m s , on the other hand, is surprisingly much better than that for m. Since the parameters of the 
clover fit are poorly determined, we do not consider them any further. 
If, instead, we fit the data keeping points at (3 > 6.0, then 



m s (M ) = 105(21) MeV[l + 1.4(1.0) GeV a] x 2 /dof = 0.3 Wilson 
m s (M^) = 104(5) MeV[l — (0.74(14) GeV a) 2 ] x 2 /dof = 1.8 Staggered 1 



For our best estimate of m s (M r j ) ) in the quenched theory we average the Wilson and staggered values 
given in Eqs. W and |TT| 



m s (MS,n = 2 GeV, Afy) = 100(21) MeV (quenched). (18) 
where we quote the largest of the extrapolation errors. 
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The rif = 2 data for m s (M^) are shown in Fig. |7|. In light of the discussion presented in 
Section |5| and |8] we only consider data with (3 > 5.5. Furthermore, the spread between the values of 
m s obtained by different collaborations at a given (3 are large, so we do not extrapolate assuming 
the expected 0(a) corrections. Instead we simply take the average in both cases. This gives 
m s = 76(10) MeV for the Wilson formulation and m s = 59(6) MeV for the staggered. For our 
final estimate in the a = limit, we take the average of these two values, 

m s (~MS, // = 2 GeV, M+) = 68(12) MeV (n f = 2) , (19) 

where the spread in the data is taken to be representative of the uncertainty in the extrapolation 
to the continuum limit. 

The estimates for m s found using Mk* match with m s (M^) as exemplified in Fig. |[ To sum- 
marize, we find that the quenched estimate of m s depends on the hadronic state (Mk versus M<j> 
or Mk*) used to extract it. While the first estimate is constrained by the lowest order chiral 
perturbation theory, the second is an independent estimate. Thus, for our final estimates we take 

m s (MS, f A = 2GeV,M 4> ) = 100 ±21 ±10 MeV (quenched), (20) 
m s (MS,n = 2 GeV,M^) = 68 ± 12 ± 7 MeV (n f = 2) . (21) 

where the second error is the 10% uncertainty in setting the lattice scale as discussed in Section §. 
Using these, one can again make a linear extrapolation in n./ to get m s ~ 55 MeV for the three 
flavor theory. However, we feel that it is important to stress that the unquenched simulations are 
still in very early stages with respect to both statistical and systematic errors. Thus, unquenched 
estimates should be considered preliminary. 

To summarize, the data show two consistent patterns. First, for a given value of a the nf = 2 
results are smaller than those in the quenched approximation. Second, the ratio m/m s (M ( j } ) is 
in good agreement with the predictions of chiral perturbation theory for both the n/ = and 2 
estimates. Quantitatively, our estimates are low and will be refined as data at more values of (3, 
possibly with n/ = 2 — 4 to bracket the physical case of n/ = 3, become available. 



7 STUDY OF THE CLOVER IMPROVEMENT 

At (3 = 6.0 there now exist data at a number of values of the clover coefficient as listed in Table ^. 
This allows us to study the effect of improvement on quark masses as a function of csw in the 
range 0-3.0. Note that the value predicted by the Alpha Collaboration for 0(a) non-perturbative 
improvement is 1.769 [11]. The data, presented in Table are plotted in Fig. H along with the 



Wilson and staggered fits reproduced from Fig. [I], and the staggered points at /3 = 6.0. Qualitatively, 
increasing the clover coefficient increases a and decreases m. In particular, the scale at csw = 1-769 
matches that from the staggered data. However, the value of m is significantly (~ 50%) different. 
For a complete 0(a) improvement with clover actions, the lattice quark mass needs an 0(ma) 
correction, i.e. m t ™P roved = m q (l ± b m m q a) where b m = —0.5 at tree level jTTJ. This is, however, 



a few percent effect only for the parameter values we have analyzed. The other two possible 
explanations are the failure of perturbation theory used to determine Z m , and possibly large 0(a 2 ) 
effects. More work, like determining Z m non-perturbatively, is required to quantify these effects. 
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Table 6: Lattice parameters of the runs at /3 = 6.0 for different values of the clover coefficient in 
the Sheikholeslami-Wohlert action. 



Label 



Ref. 



csw 



Lattice 



size 



#of 
Conf. 



Quark masses used in the fits 



Scale 1/a 
(GeV) 



a- 



[1/a) (2 GeV) 



W8 

CI 

C5 

C6 

C7 

C8 

C9 



[0 
[0 
[0 
[0 
[0 
[0 





1.0 
1.4785 
1.769 

1.92 

2.25 

3.0 



32 3 x 64 
18 3 x 64 
32 3 x 64 
24 3 x 64 
16 3 x 32 
16 3 x 32 
16 3 x 32 



170 
200 
100 
100 
100 
100 
100 



0.155, 0.1558, 0.1563 
0.1425, 0.1432, 0.1440 
0.13808, 0.13851, 0.13878 
0.1342, 0.1346, 0.1348 
0.1290, 0.1300, 0.1310, 0.1320 
0.1260, 0.1265, 0.1270, 0.1277 
0.1160, 0.1165, 0.1170, 0.1173 



2.338(43) 
1.867(81) 
2.028(96) 
1.762(78) 
1.734(26) 
1.695(42) 
1.594(48) 



0.192 
0.192 
0.192 
0.192 
0.192 
0.192 
0.192 



0.205 
0.187 
0.193 
0.184 
0.182 
0.181 
0.177 



Table 7: The results for the quenched simulations using the Sheikholeslami-Wohlert action as listed 
in Table ^. The notation is the same as in Table |j. 



Label 


k c B^ A p Bp 


m 

ma x 10 5 m 


m s (M0) 
m s a x 10 3 m s 


W8 

CI 

C5 

C6 

C7 

C8 

C9 


0.157143(03) 1.993(26) -6.663(521) 2.196(162) 
0.145483(03) 2.674(29) -6.265(1006) 1.942(288) 
0.139285(01) 2.837(50) -9.913(2491)2.866(689) 
0.135222(05) 3.268(126) -7.318(2973) 2.096(799) 
0.133062(05) 3.483(20) -8.914(243) 2.489(63) 
0.128536(07) 3.640(50) -9.934(673) 2.669(170) 
0.118151(104) 3.906(195) -13.101(714) 3.209(166) 


172(06) 5.20(11) 
201(18) 5.33(24) 
161(15) 4.89(24) 
185(17) 4.87(27) 
179(05) 4.67(07) 
179(09) 4.64(12) 
189(10) 4.67(19) 


50(05) 152(11) 
71(13) 188(27) 
45(12) 136(32) 
70(29) 183(68) 
60(02) 156(4) 
57(05) 147(9) 
51(04) 125(6) 
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8 ANALYSIS OF SYSTEMATIC ERRORS 

We have checked the validity of our analysis ignoring correlations between data by comparing 
results presented here against those obtained from a detailed error analysis using the quenched 
Wilson data obtained on 170 32 3 x 64 lattices at (3 = 6.0 P2| . For example, in Ref. [p2"| , we quote 
m B (M < f ) ) = 154(8) MeV, whereas the current analysis gives 152(11), the number plotted in Fig. [|. 
A similar agreement is seen if we compare our estimates with the results of the original analysis 
by collaborations whose data we have used as discussed in Section Thus, we believe that the 
analysis presented here gives central values in the individual points reliable to within a few percent, 
and the estimate of statistical errors to within a factor of two. 

There is accumulating evidence in high statistics data that fits linear in the quark mass, Eq. [11], 
are not sufficient, i.e. chiral logs and terms of higher order in m need to be included || [JO, H3|. 
(Chiral perturbation theory for pseudoscalar mesons does not indicate any large non-linearities as 
the change in m s /m and m u /md between the lowest and next order analyses is insignificant H.) In 



a state-of-the-art quenched calculation presented in Ref. p2| , we found that including higher order 
terms in the fits for pseudoscalar and vector mesons changes the estimates of quark masses by less 
than a few percent. In the present study we have not carried out this detailed analysis because in 
most of the data the quality of statistical precision and the number of values of light quark masses 
simulated are insufficient to uncover these effects. The deviations from linearity are found to be 
smaller for pseudoscalar mesons, which is another reason why we used them to fix m, and vector 
mesons to fix m s . In short, even though present lattice data does not rule out the possibility of 
large systematic errors due to the neglected higher order terms, there is no clear indication that 
this is the case. This issue can be addressed only when the range of quark masses used in the fits 
is increased. 

For Wilson and clover type of fermions, the errors in the calculation of k c and the k corresponding 
to light quarks is highly correlated. As a result the estimate of errors in 1/2k— 1/2k c is much smaller, 
and we find that the low statistics points give estimates that are consistent with the high statistics 
large lattice data. Thus we believe that the estimate, in the a = limit, will not change significantly 
on improving just the statistical quality of the current data. What is more important is to increase 
the range of quark masses in the chiral fits and to verify whether the lowest order correction formulae 
used in the extrapolation to a = are sufficient. We have partially addressed the second issue by 
comparing the continuum limit of data from different lattice actions. We plan to test these issues 
in the future as data at smaller m q and a become available. 

The clover action with csw = 1-0 has discretization errors starting at 0(a s a) and it is not known 
whether 0(a s a) or 0(a 2 ) errors dominate for (3 > 6.0. Nevertheless, one would have expected the 
corrections to be smaller than those for Wilson fermions. However, the current data show a near 
agreement between Wilson and clover data. With data at just three values of coupling it is hard to 
determine what form to use to extrapolate the clover action data to the continuum limit. The most 
we can say is that using a linear extrapolation, a theoretically allowed form, gives results consistent 
with the rest of our analysis. 

The lattice scale a: The value of the lattice scale a enters into the calculation in two ways. 
It is used to convert the dimensionless lattice quantity ma into physical units. A much weaker 
dependence arises when matching the lattice and continuum theories. This is because one specifies 
the renormalization point /i, at which the lattice and continuum theories are matched, in terms of 
1/a. Present lattice data show a considerable spread in the determination of a depending on the 
physical quantity used. This variation is due to a combination of discretization and quenching errors. 



For example, with our high statistics quenched Wilson data p2 |, we have calculated the lattice 
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scale fixing M p , f n , and M N to their physical values. The results are l/a(M p ) = 2.330(41) GeV, 
l/a{U) = 2.265(57) GeV, and l/a{M N ) = 2.062(56) GeV Also, the fluctuations between the 
determination of scale at a given (3 by the different groups (see Tables @, [| ||) are of the same 
magnitude. Thus, using M p to set the scale could lead to an additional 10% uncertainty in our 
estimates of quark masses, which we quote as a separate error in our final estimates. 
The rif = 2 Wilson data: We show the data for (3 > 5.3 obtained on 16 3 x 32 lattices in 
Fig. H We highlight the highest statistics points by the symbol octagon. The striking feature 
of this data is the surprisingly low value m « 1.1 MeV at (3 = 5.3. Coincidentally, the SCRI 
collaboration has reported an anomalous behavior in the calculation of the non-perturbative (3- 
function at (3 — 5.3 [j36fl . The important question is whether this is an anomalous measurement 
or is indicative of a serious problem. If this feature is borne out by future calculations, then one 
possible interpretation is that the lattice theory has a phase transition close to (3 = 5.3, similar to 
that seen in the fundamental-adjoint coupling plane of the pure gauge model. Such singular points 
are lattice artifacts, and to avoid their influence in the extrapolation to a = one has to consider 
data at weaker coupling or improve the action to stay away from such singularities. In the former 
case one would need to consider data at couplings weaker than (3 = 5.3. Motivated by the need to 
minimize contamination from such effects, we have chosen to use data at (3 > 5.5 only. 

The above scenario would also explain why the Wilson data for (3 < 5.5 (a > 0.42 GeV -1 ) appear 
to lie marginally below the staggered data, in contrast to the quenched estimates. Furthermore, 
one might worry about the presence of such artifacts at even weaker coupling. In the long run such 
a possibility is addressed by our strategy of demanding agreement, in the a = limit, of results 
obtained using different lattice discretizations of the Dirac action since the presence and influence 
of such artifacts depend on the lattice action. 



9 COMPARISON WITH PREVIOUS ESTIMATES 



The ROME/APE collaboration quotes the value m^M^) = 128(18) MeV [12| from their quenched 
simulations. Their estimate is based on an analysis of data obtained using both the Wilson and the 
clover actions at f3 = 6.0 and 6.2, but does not include a linear extrapolation to a = 0. Their number 
is consistent with the result 128(4) MeV reported by us in Ref. ||22j| . However, our present analysis 
of the global data show the presence of large 0(a) corrections for the Wilson theory. Including this 
correction results in a much lower number, m s (M^) = 88(10) MeV, after extrapolation to a = 
as shown in Fig. |j. 

The APE collaboration has recently updated their estimates using data with Wilson and Clover 
actions at f3 = 6.0, 6.2, 6.4 [pT]] . These data are included in our analysis and we agree with the raw 



numbers. However, they quote the value m s (MS,2 GeV, Mr-) = 122(20) MeV, based on taking 
an average of the data rather than doing an extrapolation in a. Since they only analyze their own 
data, they are not able to make a case for linear extrapolation as advocated by us. 



Gough et al. have also presented results for both the Wilson and clover actions for (3 



5.7, 5.9, 6.1. Even though their values for m(MS,2 GeV) at a given (3 are systematically lower 
than what we find due to the differences between the two analyses in determining a, their estimates 
of the continuum values are consistent with the results presented here. 

The analysis by Lee of the rif = 2 staggered data generated by the Columbia collaboration at 
(3 = 5.7 yielded m = 2.7 MeV |45| . This is consistent with 2.62(9), the result of our analysis of their 
data, and with the value, 2.95(14) MeV, we get for the data by Fukugita et al. |38| at the same 
coupling. 
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To conclude, the values of quark masses that we have extracted at a given lattice scale are 
consistent with previous analysis. Both the central value and the error estimates based on our 
simple analysis are consistent with the previous analysis once the differences in the regularization 
schemes are taken into account. The main new feature of our analysis is to show that the data 
roughly follow the expected 0(a) corrections for quenched Wilson and clover fermions and 0(a 2 ) 
for staggered. We also show that including these corrections give results that are consistent between 
the various discretization schemes. 



10 IMPACT ON e'/e 



The Standard Model (SM) prediction of e'/e can be written as [46 



e'/e = A jc + [c 6 5 6 1/2 + c 8 S 8 3/2 ]M r | , (22) 

where M r = (158 MeV / (m s +md)) 2 and all quantities are to be evaluated at the scale m c = 1.3 GeV. 
This form highlights the dependence on the light quark masses and the B parameters. For central 
values of the SM parameters quoted by Buras et al. j|(J, we estimate A = 1.29 x 10~ 4 , c = —1.4, 
c 6 = 7.9, c 8 = —4.0. Thus, to a good approximation e'/e oc M r ; and increases as B^ 2 decreases. 
Conventional analysis, with m s + md = 158 MeV and B X J 2 = B^ 2 = 1, gives e'/e m 3.2 x 10 -4 . 
On the other hand taking m s + pa 85 MeV, our rif = 2 estimates scaled to m c , and B^ 2 = 0.8 
2g gives e'/e « 19 x 10" 4 . This estimate lies in between the Fermilab E731 (7.4(5.9) x 10~ 4 ) g7 



and CERN NA31 (23(7) x 10 4 ) [48] measurements. Since the new generation of experiments will 



reduce the uncertainty to 1 x 10 , tests of the enhanced value are imminent. 



11 Continuum limit of the Quark Condensate 

The observation that the present lattice data for pseudoscalar mesons is well described by M% = 
B n m q allows us to calculate the chiral condensate using the Gell-Mann-Oakes-Renner relation P, [B| 

<^> GMOR = lim -ll^k. (23) 

Since m q {ipip) is renormalization group invariant, we analyze the slope of M 2 versus the renormalized 
mass at a fixed scale, i.e. m q (MS,2 GeV), for the data sets described in Tables |2] and |3|. The 
slope, after extrapolation to a = 0, directly give an estimate of 4 (ipip) / f 2 in MS scheme at ji = 



2 GeV. Our results for the quenched and dynamical configurations are displayed in Figs. |10| and |TT 
respectively. The fits, assuming a linear behavior in a for Wilson and clover and a 2 for staggered 
fermion formulations, to the quenched data with f3 > 6.0 give 



f 2 _ 

-4(#> 

ft 
-4(#> 
ft 



= 5.41(12) GeV[l- 0.72(7) GeV a] X 2 /dof = 2.6 (Wilson n f = 0) (24) 
= 5.64(15) GeV[l - 0.73(6) GeV a] x 2 /dof = Yl.\(dover n f = 0) (25) 
6.13(3) GeV[l - (0.46(2) GeV a) 2 } x 2 /dof = 19. 5 (staggered n f = 0) (26) 
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The fits for clover and staggered data do not work. In particular, the staggered data show a break 
at (3 ~ 6.0. At weaker coupling the data lie in the band 5.9(2) showing no clear a dependence. 
Taking this to be the best estimate for staggered fermions, we find that the three values are in 
rough agreement. Thus, for our final value we take the mean of these, i.e. 4(^0)//^ = 5.7(4) GeV, 
where the error covers the spread. Using the experimental value of f n = 131MeV, we then obtain 

(^0) = -0.024 ± 0.002 ± 0.002 GeV 3 (27) 
m s = -0.0023 ± 0.0006 ± 0.0004 GeV 4 (28) 

as our estimate for the quenched theory. The second error arises from the 10% scale uncertainty 
discussed in section ||. 

The behavior of the dynamical data is again not clear. Therefore, we take the average 7.3(1.0) 
of the values at the weakest coupling as our best estimate. This gives 

(foj,) = -0.031 ± 0.004 GeV 3 (29) 
m,(^) = -0.0021 ± 0.0006 GeV 4 . (30) 

The corresponding phenomenological estimates are 

= -0.0114 GeV 3 (31) 
m s $ip) = -0.0017 GeV 4 . (32) 

Thus, lattice estimates give low values for the quark masses and correspondingly high values for 
the condensate, while roughly preserving m(ipip). 



12 CONCLUSIONS: 

We have presented an analysis of m and m s using data generated by us over the years and also 
by other collaborations. The values of quark masses we have extracted for a given set of lattice 
parameters are consistent with previous analyses. Both the central value and the error estimates we 
get by reanalyzing the data are consistent with the previous reported results once the differences in 
the regularization schemes are taken into account. The main new feature of our work, based on the 
global analysis, is to show that the data roughly follow the expected 0(a) corrections for Wilson 
fermions and 0(a 2 ) for staggered. The errors in the clover action are expected to be O (a /log (a)), 
however over the limited range of a it is not surprising that the present data is consistent with a 
linear fit. What is surprising is that the 0(a) corrections are still as large as for Wilson fermions. 
The bottom line is that including these corrections give results that are consistent between the three 
discretization schemes after extrapolation to the continuum limit. 

We quote our final results in the MS scheme evaluated at /x = 2 GeV. The lattice perturba- 
tion theory is reorganized using the Lepage-Mackenzie scheme. Our best estimate of the isospin 
symmetric mass m is 3.4 ± 0.4 ± 0.3 MeV for the quenched theory. For the rif = 2 flavors there 
does not exist enough data to extrapolate to the continuum limit. The mean of the data at weakest 
couplings gives m(MS, 2 GeV) = 2.7 ± 0.3 ± 0.3 MeV. Using a linear extrapolation in rif would 
give m ~ 2.4 MeV for the physical case. 

To extract the value of m s we use the physical value of (or equivalently M^*)- Using Mk 
instead constrains m s (M^) = 25.9m since we use a linear fit to the pseudoscalar mass data. Our 
best estimate for the quenched theory is m s (MS,fi = 2 GeV, M^) = 100 ± 21 ± 10 MeV, and 



17 



68 ± 12 ± 7 MeV for the two flavor case. The variation with n/ would again suggest an even 
smaller value for the physical rif = 3 theory, however, we caution the reader that the unquenched 
simulations are still in infancy. 

In short, taking into account the various systematic errors, like allowing for a 10% uncertainty in 
the estimates due to the uncertainty in setting the lattice scale a, the present lattice results for both 
m and m s are surprisingly low compared to the numbers used in phenomenology, i.e. m = 6 — 7 MeV 
and m s = 150 — 175 MeV. The main uncertainty in the lattice results arises due to lack of control 
over the extrapolation of the two flavor data to the continuum limit, and consequently the final 
extrapolation to n/ = 3. To address these issues requires unquenched data at more values of j3. 
The data suggest that, with respect to statistical and discretization errors, the better approach is 
to use staggered fermions, however one has to confront the issue of a large Z m . On the other hand 
one needs to understand why the non-perturbatively improved Sheikholeslami-Wohlert action has 
large discretization errors. Clearly, the reliablity of these estimates will be improved in the next 
few years as more data become available. 

From a study of the variation of the pseudoscalar data as a function of the quark mass we also 
extract the quantity (ipip) / fl using the Gell-Mann-Oakes-Renner relation. After taking the contin- 
uum limit we find that the chiral condensate is roughly a factor of two larger than phenomenological 
estimates. Consequently, the estimate of the renormalization group invariant quantity m(tpip) is 
preserved. 
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Figure 1: The behavior of fn(MS,2 GeV) extracted using the quenched M n data in the TAD1 
lattice scheme defined in the text. The scale is set by M p . We do not show the fit to the clover 
data for clarity. 
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Figure 2: Three different fits to quenched Wilson data for m(MS, 2 GeV). The solid line is a linear 
fit to points at f3 > 5.93, while the dotted line includes a quadratic correction. The dashed line is 
the quadratic fit to points at f3 > 6.0. 
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Figure 3: The behavior of m(MS,2 GeV), extracted using M w data for rif — 2 simulations. The 
scale is set by M p , and the lattice scheme is TAD1. 
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Figure 4: Comparison of m s (MS,2 GeV) extracted using against the lowest order chiral pre- 
diction m s = 25.9m. The data are for the quenched Wilson theory. 
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Figure 5: Comparison of m s (MS,2 GeV) extracted using for the quenched Wilson, and stag- 
gered theories. 
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Figure 6: Results for m s (MS,2 GeV) extracted using M K , M K *, and for the quenched Clover 
action with tree-level value csw = 1 for the clover coefficient. The results for m s {M ( j > ) and m s (MK*) 
agree. We show linear fits to the rris^M^) and m s (Mjf) data. 
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Figure 7: Comparison of m s (MS,2 GeV) extracted using for the rif — 2 Wilson and stag- 
gered actions. The lines show the average value and the uncertainty in the estimate for the two 
formulations. 
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Figure 8: Behavior of m(MS, 2 GeV) versus scale a as a function of the clover coefficient csw- We 
also reproduce the fits to the Wilson and staggered data, and the staggered result at (3 — 6.0 from 
Fig. |1] for comparison. 
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Figure 9: The rif — 2 data for Wilson fermions. We show the highest statistics points at /3 = 
5.3, 5.5, 5.6 with the symbol octagon. The remaining points are shown by the symbol square. We 
also plot our estimate, the average value for (5 > 5.5, and show the errors in a for each data point. 
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Figure 10: The slope of Ml versus m q (MS, 2 GeV) for the quenched data as a function of the lattice 
spacing a(M p ). For clarity, the fit to the clover data is not shown. 
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Figure 11: The slope of M% versus m q (MS ,2 GeV) for the dynamical data as a function of the 
lattice spacing a(M p ). 
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